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1. Introduction 



At TeV colliders, the electroweak corrections are strongly enhanced by logarithms of the form 
a'ln^ with M = Mw ~ Mz and j < 21. These logarithmic corrections affect every re- 

action that involves electroweakly interacting particles and is characterized by scattering energies 
^ M^. For 2 ~ 1 TeV their impact can reach tens of per cent at one loop and several per cent 
at two loops (see for instance Refs. [|l|, ^ and references therein). Their origin is twofold. The 
renormalization of UV singularities at the scale /Xr < M gives rise to terms of the form ln(2^//i^). 
In addition, the interactions of the initial- and final-state particles with soft and/or coUinear gauge 
bosons give rise to mass-singular logarithms, i.e. ln(2^/M^) terms that are formally singular in 
the limit M — > 0. Since they originate from UV, soft and collinear singularities, electroweak log- 
arithmic corrections have universal properties that can be studied in a process-independent way 
and exhibit interesting analogies with QCD. At one loop, the leading logarithms (LLs) and next-to- 
leading logarithms (NLLs) factorize and are described by a general formula that applies to arbitrary 
Standard-Model processes [^, 0]. The properties of electroweak logarithmic corrections beyond 
one loop are investigated with two complementary approaches in the literature: (i) Evolution equa- 
tions, which are well known in QED and QCD, are applied to the electroweak theory in order to 
obtain the higher-order terms through a resummation [Q, ^, ^, ||, ^]. In this approach the evolution 
is split into two regimes, where the electroweak interaction is described in terms of SU(2)xU(l) 
and Uem(l) symmetric gauge theories.' This splitting is assumed to correctly reproduce all relevant 
implications of electroweak symmetry breaking in LL, NLL and NNLL approximation, (ii) This 
assumption and the resulting resummation prescriptions can be checked and improved by means 
of explicit diagrammatic two-loop calculations based on the (spontaneously broken) electroweak 
Lagrangian. So far, all existing diagrammatic results are in agreement with the resummation pre- 
scriptions, however up to now only a small subset of logarithms and processes has been computed 
explicitly. While the LLs [|lj] and the angular-dependent subset of the NLLs [|l^] have been de- 
rived for arbitrary processes, complete diagrammatic calculations at (or beyond) the NLL level 



exist only for matrix elements involving massless external fermions [13, 14, 15]. In the literature, 
no explicit NLL calculation exists for reactions involving massive scattering particles. 

In the following we present a recent two-loop NLL calculation for n-fermion processes 
/1/2 h - --fn involving an arbitrary number of leptons and/or quarks. While all fermions with 
nif <^M are treated as massless, the top-quark mass effects are taken into account. For all invariants 
= {Pi + Pj)^ we assume |r,j | ~ 2^ ^ M^, and the differences between the W, Z, H and t mass are 
taken into account. Soft and collinear singularities from massless virtual photons are regularized 
dimensionally and arise as £-poles in D = 4 — 2£ dimensions. For consistency, the same power 
counting is applied to ln(2^/M^) and l/e singularities. Thus, in NLL approximation we include 
all £^'^W^*^(2^/M^) terms with total power j = 2,1 at one loop and 7 = 4, 3 at two loops. The 
photonic singularities are factorized in a gauge-invariant electromagnetic term. The remaining part 
of the corrections - which is finite, gauge invariant, and does not depend on the scheme adopted to 
regularize photonic singularities - contains only ln{Q^/M^) terms. The divergences contained in 
the electromagnetic term cancel if real-photon emission is included. 



'a new approach based on soft-collinear effective theory has been developed in Ref. |p^. 
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2. Extraction of soft, coUinear and UV contributions from Feynman diagrams 



Let us consider the two-loop diagrams that produce NLLs within the 't Hooft-Feynman gauge. 
These diagrams are obtained from the LL one-loop diagrams, i.e. diagrams with external-leg ex- 
change of soft-coUinear gauge bosons, via insertions of one-loop sub-diagrams that produce an 
additional logarithm. We start with the diagrams that do not contain Yukawa interactions. 




^2 





■ (2.1) 



Here V,- = A,Z,VK*, and the external lines i,j,k = l,...,naie massless and/or massive fermions. 



In the first four diagrams in (gjj) the gauge boson V2 can couple to an external line or to an internal 
propagator. In the latter case the only relevant region is the one where V2 is collinear to an external 
fermion and, using collinear Ward identities, we find that this type of contributions cancel [[I^, p^ ]. 
The remaining diagrams contain a tree sub-diagram that corresponds to the original process and 
depends only on the external momenta.^ These diagrams are called factorizable, since their NLL 
contribution factorizes into the w-fermion tree amplitude times a two-loop NLL factor. For instance. 



the factorizable diagrams associated with the fourth diagram in (2.1) are 






(2.2) 



To extract the NLLs of soft/coUinear origin, we exploit the fact that in the soft/coUinear limit the 
coupling of a gauge boson V to the i-th external fermion yields a factor —2eiy{pi + q)^, where 
q and pt are the gauge-boson and fermion momenta, while if is the gauge-group generator. For 
instance, for the last diagram in ( |2.2[ ) we obtain 



4e 



2£ 



-^2)^:3+ perm.] /f'/f/g^^ 



(2;r)2^ {q\ 



-Ml){q 



Mi){qi-M^^){lf-mi){q-mi)m 



where q^ = —qi —q2,h = Pi — q\,h = Pk — qi and Z3 = pj — q^. The «-fermion tree amplitude ^0 
must be regarded as a vector carrying SU(2) indices associated with the isospin of the scattering 
fermions, and the generators if act as SU(2) matrices on the i-th fermion. 

In addition to soft/collinear contributions, there are also two-loop NLL terms that arise from 
UV divergences in the one-loop sub-diagrams. In this case, the UV singularities are removed by 
means of a minimal subtraction at the scale jj.^ = Q^. As a consequence, only those UV-divergent 
sub-diagrams characterized by scales of order <C produce NLL terms. In particular, non- 
factorizable diagrams do not produce NLL terms of UV origin. The explicit factorization of the 
NLL contributions of UV origin is obtained by means of projector techniques [15, |^ since the 



^In these tree sub-diagrams, which are denoted with an 'F' in (2.2), the loop momenta must be set to zero. This 
follows from the collinear Ward identities. 
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soft-coUinear approximation is not applicable in this case. After the abovementioned minimal UV 
subtraction, we perform a finite renormalization that restores the on-shell normalization of the wave 
functions and translates the couplings to the minimal-subtraction scheme at the scale /Xr. 

For processes involving heavy quarks, in addition to the contributions listed above also dia- 
grams with Yukawa interactions contribute. Many of them are suppressed due to the behaviour of 
the interactions of scalar fields in the soft/coUinear and M —fO hmits. Only three types of diagrams 
are not suppressed. 






(2.3) 



Moreover it turns out that the sum of the NLL terms resulting from these diagrams cancels owing 



to global gauge invariance of the Yukawa sector ||160. As a consequence, Yukawa interactions 
contribute only through fermionic self-energy diagrams that appear in the renormalization of the 
heavy-quark wave functions. 



3. High-energy expansion of the loop integrals 

The loop integrals have been calculated in the » limit to NLL accuracy using the 
sector-decomposition technique JT7| , p^ ] and, alternatively, the method of expansion by regions 
[1191,0,0,01. 



3.1 Sector decomposition 

Sector decomposition is based on the Feynman-parameter representation. Exploiting general 
properties of soft, collinear and UV singularities, one can factorize these singularities in Feynman- 
parameter space in such a way that the divergent integrations assume a standard form. If all parti- 
cles are massless these singular integrals are trivial and produce simple £-poles [[l7|]. In presence 
of massive particles, the general solution is available to NLL accuracy JT^]. This permits to con- 
struct, in a completely automatized way, finite integral representations for the coefficients of all 
NLL terms of type e^'^ln-'^'^^Q^ /M^). The integrands are smooth and can easily be integrated nu- 
merically.^ Moreover, although these integrals do not have a standard form, it turns out that they 
are relatively simple. In practice they can be solved in analytic form by means of standard tricks 
(partial fractioning, integration by parts, etc.) combined with a table of elementary integrals. Also 
this part of the calculation is highly automatized, however, the set of integration rules might need 
to be extended when new diagrams are computed. 

3.2 Expansion by regions 



Within the expansion-by-regions method [19, 20, 21, E3], the integration domain of the loop 



momenta is divided into regions corresponding to the asymptotic limit considered. The integrand is 
appropriately expanded in every relevant region, and each of the expanded terms is integrated over 



This holds only in the unphysical region where all invariants r,y are negative, but allows for very useful checks. 
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the whole integration domain. Each expanded term has a unique order in powers of Q and M, but the 
on-shell momenta pi of massive fermions involve two scales, pj ~ and Ipipj --^ Q^. In order to 
separate these scales, the massive momenta are reparametrized in terms of light-like momenta pij 
as Pi = Pi + p} / {2piPj)pj. For the loop integrals needed in this calculation, the following regions 
for each loop momentum k are relevant: hard, soft and ultrasoft regions, where is of the order Q^, 
and M^/ Q^, respectively; coUinear, softcoUinear and ultracoUinear regions, where k is coUinear 
to one of the external momenta and k^ is of the order M^, M^/Q^ and M^/Q^, respectively. The 
latter two regions contribute only if the external fermions are massive. The expanded loop integrals 
have been evaluated using Mellin-Bames representations, from which the extraction of NLLs has 
been automatized. In addition, Mellin-Barnes representations of the unexpanded integrals have 
been used to check the completeness of the set of regions used for the expansion. 



4. Results and conclusions 

The renormalized 2-loop NLL amplitude is expressed in terms of combinations of SU(2) ma- 
trices acting on the tree-level amplitude. This expression is simplified by means of SU(2) algebraic 
identities and additional relations resulting from global gauge invariance. The final result is con- 
sistent with the double-exponentiation formula^ 



■^ = expl^I^lexp(^i:5rl(l + ^I^ 



8^ 



(4.1) 



The term 5?™ behaves as if all electroweak gauge bosons would have the same mass M and is 
symmetric with respect to SU(2)xU(l) transformations, 

- £ /f/J/(£,M,-r,;)--^ 



sjsew 

"7 



(47r) 



8\ 



2^'^; > (1) , „27.fl7.a, (1) 



b\'>+giTf'Tfb'^>\j{eM,^ii)- (4.2) 



Here T" and Yi are the generators of the SU(2)xU(l) group, gi 2 are the corresponding coupling 
constants, and b^p2 the one-loop j8 -function coefficients [15|. The NLL functions read 



2 1 

I{e,M,r) = -L}--L^e--L^e^ + 



3 r\ yf g\m\ 

.2 Vev crse^Mi, 



and J{e,M,^l) = [l{2e,M,Q^) - {Qy^lyi{e,M,Q^)] /e with L = IniQ^/M^). In Q q 



Y^v^i u ' ^^'^ the Yukawa terms proportional to mf contribute with factors y]^^ = l,y^ =2,y^ = 0, 
for heavy quarks / = b,t with chirality K",- = R,L, and yf = otherwise. The term dfj^ in ( plj ) 
depends on the 7-W mass gap and contains all soft/collinear singularities due to massless pho- 
tons. These singularities involve single and double £-poles. Thus the terms df-^^ and dfj in ( pTT| ) 
need to be expanded up to ^(e^). We find that d^p behaves as in QED [[l6|]. Finally we find that 
aU contributions depending on the Z-W mass difference can be factorized into the one-loop term 
= -if if In (Ml/Mlf) [2L + 2L^e + L^e^] . Apart from these latter ln(M|/M^) terms, all effects 
of symmetry breaking (terms proportional to Mw/Mz, mixing parameters and vacuum expectation 



*In the following we set fi^ = Q^e^ / (Atz) 
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value) disappear due to non-trivial cancellations between different Feynman diagrams. These can- 
cellations originate from relations between the gauge-boson masses, the mixing parameters, and 
the vacuum expectation value. The result is in agreement with the resummation prescriptions pro- 
posed for massless |8]] and massive fermions [Q] and exhibits strict analogies with the general 
form of two-loop singularities in massless QCD [^]. 
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